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SUMMARY 


The  state  of  the  art  for  assessment  of  defect  criticality  in 
structural  composite  laminates  has  been  reviewed  and  summarized. 
The  significance  of  related  ongoing  programs  has  been  reviewed 
through  a  literature  survey  and  the  organization  of  a  limited 
attendance  symposium  on  the  subject. 

Tests  have  been  performed  to  determine  the  range  of  validity 
of  criticality  criteria  for  disbonds  in  laminated  beam  and  plate 
type  structures,  which  were  developed  in  previous  related  programs 
Data  correlation  studies  have  shown  the  usefulness  of  a  linear 
elastic  fracture  mechanics  approach,  of  methods  of  stress  analy¬ 
sis  based  on  2-rD  elasticity  and  modified  laminated  plate  theories, 
as  well  as  of  semi-empirical  growth  laws  for  cyclic  loading,  for 
assessing  growth  of  disbonds  under  transverse  shear. 
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FOREWORD 


This  report  summarizes  the  work  done  for  the  Naval  Air 
Development  Center,  Warminster,  PA  18974,  under  contract 
N  62269-80-C-0271  during  the  period  September  30,  1980  to 
August  5,  1981. 
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thanks  are  due  to  Mr.  Robert  A.  Blake,  Jr.  for  his  extensive 
work  in  the  experimental  studies. 
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Analytical  Evaluation  of  Criticality  of  Defects  in  Structural 
Composite  Laminates,"  which  was  sponsored  under  the  contract 
and  was  held  in  Philadelphia  on  November  24-25,  1980,  is  a  part 
of  the  final  report  documentation.  These  proceedings  are  con¬ 
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INTRODUCTION 


Nondestructive  evaluation  of  criticality  of  defects  in  a  com¬ 
posite  structure  or  a  structural  component  poses  a  challenging  prob¬ 
lem  to  the  analyst.  The  first  step  for  such  an  evaluation  is  de¬ 
termination  of  the  size,  shape  and  location  of  an  isolated  defect 
or  of  a  number  of  closely  spaced  defects  by  nondestructive  inspec¬ 
tion  techniques.  The  next  step  is  to  assess  residual  stiffness, 
strength  or  life  in  the  presence  of  these  defects  by  the  use  of 
appropriate  analytical  methods  for  stress  analysis ,  as  well  as  to 
predict  failure  and  the  expected  growth  mechanism  of  the  defects . 

At  the  present  time  only  limited  use  is  made  of  NDE  concepts  for 
quantitative  evaluation  of  residual  properties  because  of  the  limita 
tions  of  contemporary  inspection  techniques  and  analytical  methods 
of  assessing  flaw  criticality. 

Work  under  three  previous  contracts,  N62269-77-0092 ,  N62269- 
78-C-0111  and  N62269-79-C-0209 ,  were  directed  towards  the  develop¬ 
ment  of  a  NDE  methodology  for  assessing  criticality  of  isolated  dis¬ 
bonds  in  one-  and  two-dimensional  laminated  structural  members  un¬ 
der  transverse  shear  and/or  compressive  loads.  These  investigations 
consisted  of  analytical  modeling,  as  well  as  nondestructive  and 
destructive  tests  to  establish  criticality  criteria  as  functions  of 
stress  level,  fatigue  cycles  and  environment.  Mechanisms  of  growth 
of  flaws,  and  failure  predicted  by  analytical  methodology,  were  cor¬ 
related  with  experimental  data  with  reasonable  accuracy.  Further 
tests  were,  however,  deemed  necessary  to  establish  the  range  of 
validity  of  the  proposed  methodology. 

Subsequent  to  the  start  of  work  under  the  three  investigations 
discussed  above,  other  programs  were  initiated  in  the  general  area 
of  defect  significance  in’ composite  materials.  It  appeared  worth¬ 
while  to  review  the  approaches  taken  in  such  programs  for  assess¬ 
ment  flaw  criticality.  Therefore,  the  objectives  of  the  current 
work  were: 
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1.  To  summarize  and  assess  significant  on-going  programs 
in  the  area  of  flaw  criticality  in  composites  so  as  to  obtain  an 
integrated  picture  of  the  current  state  of  the  art  in  this  field. 

2.  To  conduct  a  number  of  critial  tests  and  data  correlation 
studies,  to  establish  the  range  of  validity  of  the  criticality 
criteria  developed  during  the  previous  investigations  as  applied 

to  disbonds  with  1-D  and  2-D  planforms. 
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BACKGROUND 


In  order  to  arrive  at  a  definition  of  a  critical  flaw,  it 
is  necessary  to  formulate  an  analytical  model  which  is  capable 
of  predicting  the  direction  and  rate  of  its  growth  and  its  crit¬ 
ical  size  (at  failure) .  Due  to  the  complexity  of  the  problem 
of  failure  in  composite  laminates,  it  is  impossible  for  a  single 
analytical  model  to  encompass  the  whole  spectrum  of  fracture  and 
fatigue  behavior.  In  other  words,  different  analytical  models 
should  be  developed  to  predict  growth  of  different  types  of  flaws 
such  as  surface  flaws,  through-the-thickness  notches  and  delami¬ 
nations. 

An  initial  investigation  of  flaw  criticality  (ref.  1)  was 
concerned  with  slit  notches  and  inter laminar  disbonds  in  short 
beams  under  static  and  fatigue  loads.  The  analysis  employed  for 
interlaminar  disbonds  was  of  the  "strength  of  materials"  type. 

The  values  of  the  specific  adhesive  surface  fracture  energy,  as 
determined  from  correlation  of  test  data  for  quasi-static  growth 
of  disbonds  with  analytical  calculations,  showed  some  discrep¬ 
ancy.  For  this  reason,  a  more  sophisticated  two-dimensional 
stress  analysis  for  prediction  of  quasi-static  disbond  propaga¬ 
tion  was  later  developed  (ref.  2).  It  was  observed  that  the 
simple  "strength  of  materials"  solution  yields  estimates  of 
failure  load  close  to  those  obtained  by  the  use  of  2-D  elas¬ 
ticity  theory  for  sufficiently  large  values  of  the  ratio  of  dis¬ 
bond  length  to  laminate  thickness.  The  simple  solution,  however, 
becomes  increasingly  nonconservative  as  the  disbond  length  is  re¬ 
duced.  In  the  latter  study,  the  stress  intensity  factors  cal¬ 
culated  from  the  2-D  elasticity  solution  were  used  to  demon¬ 
strate  the  usefulness  of  the  commonly  used  semiempirical  crack 
growth  law  for  predicting  growth  of  disbonds  under  cyclic  shear. 
Test  data  from  a  series  of  tests  and  curve  fitting  procedures 
were  employed  to  obtain  the  two  constants  used  in  the  crack 
growth  law.  A  criterion  for  criticality  of  disbonds  in  a 
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laminate  under  compressive  loads  was  also  developed  and  tested 
with  the  help  of  experimental  data.  This  criterion  predicts  a 
buckling  type  of  failure  of  unsupported  parts  of  the  compression 
skin. 

The  criticality  of  disbonds  with  circular  and  nearly  ellip¬ 
tic  planforms  in  laminated  plates  under  transverse  shear  and 
compressive  loads  were  considered  in  reference  3.  As  in  the 
case  of  1-D  disbonds,  the  failure  under  compressive  loads  was 
found  to  be  due  to  buckling  of  the  unsupported  skin.  Growth  of 
disbonds  under  cyclic  shear  was  found  to  occur  in  a  slow  regu¬ 
lar  pattern.  An  analytical  methodology  based  on  a  modified 
laminated  plate  theory  was  developed  which  yields  flaw  criti¬ 
cality  criteria  and/or  growth  laws  under  compressive  loads  as 
well  as  quasi-static  and  cyclic  shear  stresses.  Because  of  the 
nature  of  the  theory  being  used  for  stress  analysis,  the  theory 
is  likely  to  yield  higher  values  of  failure  load  for  disbonds 
of  very  small  size.  Further,  in  a  majority  of  tests  conducted, 
the  flaws  were  located  in  the  mid-plane  of  the  laminate  and  be¬ 
tween  two  layers  of  the  same  orientation.  Therefore,  more  tests 
were  conducted  in  the  present  study  to  examine  the  range  of 
validity  of  the  proposed  criticality  criteria  as  regards  to  its 
location,  size,  shape  and  sharpness  of  the  disbond  tips. 

Several  other  current  studies  in  the  area  of  defect  signif¬ 
icance  (ref.  4  for  example)  have  some  aspects  in  common  with  the 
investigations  discussed  above,  although  the  dominant  stress- 
states  and  other  details  in  such  studies  are  different.  It  is, 
therefore,  of  interest  to  assess  the  significance  of  results  of 
such  programs  as  applied  to  defects  of  the  type  under  considera¬ 
tion.  Such  an  assessment  was  attempted  in  this  study  through 
literature  search  and  organization  of  a  symposium  on  defect 
criticality  in  structural  composite  laminates. 
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LITERATURE  SEARCH 


A  literature  search  on  Defect  Criticality  in  Composite  Ma¬ 
terials,  covering  a  period  of  the  past  three  years,  was  conducted 
through  the  U.S.  Small  Business  Administration  by  the  Delaware 
Technical  and  Community  College.  NTIS  and  NASA  data  bases  were 
utilized  in  this  search.  In  addition  to  this  search,  a  study 
of  recent  relevant  works  was  also  conducted.  Results  are  given 
in  the  attached  bibliographies  B1  and  B2.  Bibliography  B1  lists 
the  works  related  to  criticality  of  through-the-thickness  flaws 
(cracks,  notches  or  holes)  under  static  and  cyclic  loading.  It 
is  known  that  these  kinds  of  flaws  give  rise  to  interlaminar 
damages  as  well  as  matrix  or  fiber-matrix  interface  damages  par¬ 
allel  to  the  fibers  in  some  of  the  laminae  in  a  laminated  com¬ 
posite.  Therefore,  a  majority  of  these  studies  is  focused  on 
study  and  detection  of  such  damage  growths.  Use  of  fracture 
mechanics  methodologies  and  semiempirical  crack  growth  laws  as 
well  as  crack  growth  in  unidirectional  and  random  fiber  rein¬ 
forced  composites  are  considered  in  other  studies.  Use  of  NDI 
is  reported  in  a  few  studies  (B1.15,  1.17  for  example). 

Bibliography  B2  lists  works  related  to  delamination  growth 
in  composite  materials.  B2.1  lists  the  studies  on  delamination 
growth  in  the  Mode  I  (peel)  mechanism.  Test  «!>«thods,  as  well  as 
applications  to  static  and  cyclic  loading,  are  examined.  Growth 
of  delaminations  in  Mode  II  alone  is  considered  in  very  few 
studies,  as  listed  in  B2.2.  Test  specimen  designs  for  Mode  II 
and  Mode  III  extensions  are  also  attempted  in  B2.1.7. 

Combined  mode  effects  due  to  the  presence  of  surface  flaws, 
notch  or  free  edges,  are  considered  in  several  works  listed  in 
B2.3.  All  of  the  works  listed  in  B2.1,  B2.2  and  B2.3  are  based 
on  calculation  of  strain  energy  release  rates  and/or  stress  in¬ 
tensity  factors.  A  majority  of  the  studies  listed  in  B2.3  makes 
use  of  2-D  finite  element  calculations  and  some  of  these  utilize 
special  2-D  crack  tip  elements  near  the  delaminations.  The 
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study  B2.1.7  makes  use  of  3-D  finite  elements.  Transverse  crack¬ 
ing  in  90°  layers  is  examined  in  B2.3.1  and  B2.3.7. 

Effects  of  fatigue  and  compression  loads  are  examined  in  the 
works  listed  in  B2.4.  Debonding  between  metals  and  composites 
are  considered  in  B2.4.5,  B2.4.6  and  B2.4.8.  Other  studies  show 
that  fatigue  loading  induces  a  considerable  amount  of  delamina¬ 
tion  growth  before  final  failure  may  be  expected.  Under  compres¬ 
sion  fatigue,  final  failure  is  expected  due  to  buckling  of  dis¬ 
bonded  skins.  Calculation  of  buckling  loads  for  disbonded  com¬ 
pression  skins  is  also  reported  in  B2.2.1. 

A  considerable  amount  of  experimental  work  is  reported  in 
many  of  the  studies  for  determination  of  fracture  toughness  or 
critical  strain  energy  release  rates  and  growth  laws  under  cyclic 
loading  as  well  as  residual  strength.  Use  of  nondestructive  in¬ 
spection  techniques  is  also  reported  in  some  studies. 
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SYMPOSIUM  ON  DEFECT  CRITICALITY 


A  limited  attendance  symposium  on  "Nondestructive  and  Ana¬ 
lytical  Evaluation  of  Criticality  of  Defects  in  Structural  Com¬ 
posite  Laminates"  was  organized.  The  meeting  was  held  on  Novem¬ 
ber  24-25,  1980  in  Philadelphia.  Twenty-seven  research  workers 
active  in  the  areas  of  materials  and  structures  attended  the 
symposium.  Abstracts  and  vu-graphs  of  presentations,  as  well  as 
the  comments  and  discussions  by  the  attendees,  can  be  found  in 
the  informal  Symposium  Proceedings,  MSC  TPR  1107/1158.  The  sym¬ 
posium  was  arranged  for  the  purpose  of  reviewing  the  state  of  the 
art  of  defect  criticality  assessment  and  providing  a  forum  for 
discussion  on  future  research  needs  in  the  disciplines  of  struc¬ 
tural  analysis  as  well  as  materials  development.  Presentations 
and  discussions  dealt  with  a  variety  of  analytical  and  experimen¬ 
tal  tools  for  assessment  of  defect  criticality.  In  addition,  two 
sessions  were  devoted  to  general  discussions  on  needs  and  capa¬ 
bilities  as  well  as  recommendations  for  future  research. 

Most  of  the  presentations  focused  on  ongoing  research  ac¬ 
tivities.  Work  on  various  types  of  defects  were  reported,  and 
the  area  of  greatest  activity  appeared  to  be  the  growth  of  de¬ 
laminations  under  quasi-static  and  cyclic  loading.  Delaminations 
of  comparitively  large  size  as  well  as  small  inherent  flaws  near 
edges  and  those  inside  a  composite  were  considered.  Studies  on 
dispersed  voids  and  through-the-thickness  defects  were  also  re¬ 
ported.  The  presentations  contained  in  the  proceedings  are  sum¬ 
marized  below  for  the  purpose  of  assessing  the  state  of  the  art, 
which  is  discussed  in  the  next  section. 

1.  "Delamination  in  Graphite/Epoxy  Composites"  by  D.  J. 
Wilkins.  This  work  deals  with  the  use  of  3-D  brick 
finite  elements  for  calculation  of  strain  energy  re¬ 
lease  rates  for  delamination  type  defects  in  aircraft 
structural  components  as  well  as  test  coupons  de¬ 
signed  for  measurement  of  critical  values  of  Mode  I 
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and  Mode  II  energy  release  rates.  Growth  laws  under 
cyclic  loading  are  also  determined  experimentally. 

Use  of  such  laws  in  predicting  damage  growth  under 
fatigue  loading  is  also  considered. 

2.  "Flaw  Criticality  in  Graphite/Epoxy  Composites"  by  R. 

B.  Pipes  and  "Analytical  Modeling  of  Delaminations  in 
!  Laminated  Composites"  by  S.  N.  Chatterjee.  These  works 

deal  with  the  studies  conducted  by  Materials  Sciences 
Corporation  under  two  previous  contracts,  N62269-79-C- 
0209  and  N62269-78-C-0111.  Use  of  ultrasonic  C-scans 
I  for  monitoring  damage  growth  is  considered  in  the  first 

work.  In  the  second  study,  2-D  elasticity  theory  is 
used  for  through-the-thickness  disbonds  in  beam-type 
structures  and  a  modified  laminated  plate  theory  is 
1  utilized  for  elliptic  disbonds  in  plate  type  structures 

for  the  purpose  of  strain  energy  release  rate  calcula¬ 
tions.  Experimental-analytical  data  correlations  are 
performed  for  determining  Mode  II  critical  energy  re- 
|  lease  rate  and  crack  growth  laws  under  cyclic  loading. 

Quasi-static  buckling  loads  for  disbonded  compression 
skins  are  also  calculated  and  compared  with  experimen¬ 
tal  data. 

|  3.  "Investigation  of  Damage  Growth  in  Composites"  by  M.  M. 

Ratwani.  This  work  is  based  on  the  fact  that  inherent 
disbonds  usually  grow  under  fatigue  loading  and  it  is 
expected  that  such  disbonds  will  grow  in  the  interface 
<  where  the  transverse  shear  stress  is  a  maximum  and  ex¬ 

ceeds  a  threshold  value.  Such  shear  stresses  near  a 
hole  or  a  free  edge  are  calculated  by  a  finite  element 
method  on  the  assumption  that  the  flaws  are  absent. 
Although  the  growth  of  the  disbonds  is  usually  two- 
dimensional,  it  is  postulated  that  it  can  be  charac¬ 
terized  by  a  crack  growth  law  of  the  self-similar  type. 
Residual  strength  and/or  lifetime  under  constant  ampli- 
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tude  loading  or  spectrum  fatigue  loading  are  calcu¬ 
lated  on  the  assumption  that  buckling  is  the  mechanism 
of  final  failure.  Buckling  stress  is  computed  by  a 
formula  similar  to  the  Euler  type  column  buckling  equa¬ 
tion  with  a  multiplier  and,  therefore,  the  critical 
value  of  the  characteristic  dimension  of  the  flaw  can 
be  computed.  Some  of  the  parameters  in  the  model  are 
chosen  to  fit  experimental  results  from  compression 
fatigue  tests.  Attempts  are  made  to  monitor  damage 
growth  by  x-ray  radiography  and  C-scans. 

4.  "The  Criticality  of  Free  Edge  Delaminations  in  Graphite/ 
Epoxy  Laminates"  by  A.  S.  D.  Wang.  Strain  energy  re¬ 
lease  rates  are  calculated  by  the  use  of  2-D  finite 
elements.  Possibilities  of  delamination  growth  near 
free  edges,  as  well  as  near  transverse  cracks  in  the 
90°  layers,  are  examined  and  the  sequence  of  the  growth 
of  transverse  cracks  and  delaminations  is  postulated 
for  some  laminate  configurations. 

5.  "Flaw  Induced  Delaminations  in  Fatigued  Graphite/Epoxy 
Composites:  Effect  on  Lifetime  and  Ultrasonic  Detec¬ 
tion"  by  R.  J.  Richards.  It  is  shown  in  this  study 
that  delamination  fatigue  failures  in  laminated  com¬ 
posites  initiate  at  intrinsic  laminar  flaws.  Ultra¬ 
sonic  detection  of  flaws  is  correlated  with  optical 
observations.  An  attempt  is  made  to  find  a  relation¬ 
ship  between  lifetime  and  the  stress  intensity  factor 
range,  AK  . 

6.  "Development  of  Acceptance  Criteria  for  Graphite/Epoxy 
Structures"  by  B.  L.  Riley.  This  is  a  study  on  the 
effect  of  porosity  on  failure  under  combined  action  of 
compressive  and  shear  stresses.  Naturally  and  artifi¬ 
cially  induced  porosities  are  compared  and  a  test  pro¬ 
gram  is  carried  out  to  determine  the  failure  criteria 
under  combined  stresses. 
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"Progressive  Failure  Strength  Model  for  Advanced  Com¬ 
posite  Laminates"  by  D.  Y.  Konishi.  Progressive  lam¬ 
ina  by  lamina  failure  near  a  crack  or  a  hole  is  ex¬ 
amined.  Stresses  are  computed  by  anisotropic  elas¬ 
ticity  theory.  Failure  in  a  lamina  is  assumed  to  occur 
when  the  stresses  at  a  characteristic  dimension  away 
from  the  crack  tip  or  hole  boundary  fall  on  a  failure 
envelope.  Stiffnesses  of  the  laminate  are  adjusted 
depending  on  fiber  or  matrix  failure  in  the  failed  lam¬ 
ina.  Failure  of  the  laminate  occurs  when  the  principal 
load  bearing  lamina  fails.  Analytical  predictions  are 
compared  with  experimental  data. 

"Nondestructive  and  Analytical  Evaluation  of  Criticality 
of  Defects  in  Structural  Laminates  -  An  Airframe  Con¬ 
tractor's  Perspective"  by  J.  A.  SuareTzT"  The  approaches 
for  nondestructive  inspection  and  analytical  as  well  as 
experimental  evaluation  of  criticality,  which  are  cur¬ 
rently  being  used  or  are  being  considered  for  use  in  the 
near  future,  are  examined.  Acoustic  emission  appears  to 
be  a  useful  tool  for  detecting  and  locating  fatigue  dam¬ 
age.  The  effect  of  notches  on  static  strength  is  con¬ 
sidered  to  be  a  phenomenon,  for  which  analytical  methods 
have  been  validated.  Structural  testing  is  considered 
essential  for  assessing  criticality,  but  the  need  for 
analysis  methods  to  reduce  testing  costs  is  emphsized. 
"Investigation  of  the  Response  of  Composite  Materials 
by  Nondestructive  Evaluation"  by  E.  G.  Henneke.  Use 
of  various  NDI  techniques  are  examined  and  damage  ac¬ 
cumulation  in  composites  under  static  and  cyclic  load¬ 
ing  are  studied.*  Attempts  to  obtain  a  relationship  of 
characteristic  damage  states  to  engineering  response 
are  described. 
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Overviews  of  sponsored  research  by  the  Army,  Navy,  AFWAL 
and  NASA  describing  results  and  objectives  of  various  programs 
were  presented  by  R.  J.  Shuford,  D.  Mulville,  G.  P.  Sendeckyj 
and  G.  L.  Roderick,  respectively.  A  special  topic  of  interest 
is  the  use  of  a  large  deformation  theory  in  characterizing  de¬ 
lamination  growth  under  compression  loads  by  workers  at  NASA. 
General  discussions  indicate  that  delamination  is  the  most  cri¬ 
tical  type  of  defect  in  composite  materials,  and  in  aircraft 
structural  components  the  critical  areas  are  those  where  there 
are  sharp  changes  in  geometry  like  free  edges,  ply  drops,  etc. 
Delaminations  created  inside  a  laminate  by  impact,  like  tool 
drops,  etc.,  and  those  initiated  from  inherent  defects,  can  al¬ 
so  be  critical. 
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DISCUSSION  ON  STATE  OF  THE  ART 


The  following  is  a  discussion  on  the  current  state  of  the 
art  of  defect  criticality  assessment  as  it  appears  from  the  lit¬ 
erature  search  and  the  presentations  at  the  symposium. 

Ultrasonic  C-scans  and  Optical  as  well  as  Visual  Inspection 
appear  to  be  the  mose  widely  used  tools  for  NDI  at  this  point  in 
time.  X-ray  Radiography  and  Accoustic  Emission  have  good  poten¬ 
tials  in  this  direction.  Other  techniques  like  Thermography ,  Dye 
Penetrants  and  Attenuation  Measurements  are  still  in  the  stage  of 
research  and  development.  Very  few  of  the  NDI  techniques  have 
been  shown  to  yield  accurate  quantitative  measurements  of  defect 
size  and  shape  (see  B.1.15,  1.17  and  MSC  TPR  1107/1158).  Others 
are  useful  for  detecting  and  locating  the  flaws. 

Although  various  analytical  models  and  criteria  have  been 
developed  for  through-the-thickness  defects  like  cracks  or  holes 
and  the  damage  mechanisms  which  develop  under  static  and  cyclic 
loading  are  relatively  clear,  quantitative  assessment  of  critical¬ 
ity  of  such  defects  remains  a  complicated  task.  The  reason  for 
this  complexity  is  the  growth  of  delaminations  between  various 
laminae  and  disbonds  in  the  matrix  and/or  fiber  matrix  inter¬ 
faces  in  the  individual  layers.  For  this  reason,  use  of  frac¬ 
ture  mechanics  methodologies  to  this  problem  has  met  with  partial 
success.  Under  static  loading,  use  of  a  stress  criterion  at  a 
point  a  little  distance  away  from  the  cracktip  (ref.  5)  is  con¬ 
sidered  as  a  viable  method  by  many  investigators  as  well  as  in 
many  industries  (see  presentation  by  J.  A.  Suarez  in  symposium 
proceedings) .  Finite  element  modeling  or  progressive  damage 
growth  modeling  as  reported  by  D.  Y.  Konishi  in  the  symposium 
may  yield  more  realistic  assessments  of  damage  growth,  although 
they  have  their  own  limitations,  especially  from  the  point  of 
accuracy  as  well  as  cost.  Modeling  of  damage  growth  under  cyclic 
loading  is  still  in  the  state  of  research  and  development  (see 


NADC-80048-60 


bibliography  B-l) .  There  exists  a  strong  feeling,  however,  that 
tension-tension  fatigue  is  not  very  critical  for  fiber  composite 
laminates.  Compression  fatigue  is  critical  and  it  appears  that 
delamination  growth  between  the  laminae  can  be  modeled  by  the 
use  of  equations  similar  to  the  semiempirical  crack  growth  laws 
commonly  used  for  metals  (ref.  4).  The  usefulness  of  this  ap¬ 
proach  has  been  demonstrated,  but  a  considerable  amount  of  ground¬ 
work  is  still  to  be  performed  to  put  the  analytical  methodology 
for  assessing  two-dimensional  growth  of  disbonds  starting  from  the 
edge  of  a  hole  or  straight  boundary  on  a  sound  footing. 

Growth  of  disbonds  is  an  important  and,  in  the  opinion  of 
many  investigators,  the  most  critical  phenomenon,  not  only  in  com¬ 
pression  fatigue  but  also  in  many  other  circumstances  like  the  one 
of  transverse  shear  loading,  where  catastrophic  growth  can  occur 
under  quasi-static  loading.  Growth  of  isolated  disbonds  in  Mode 
I,  Mode  II  or  combined  modes  under  quasi-static  or  cyclic  loading 
has  been  modeled  by  various  investigators  using  fracture  mechanics 
methodologies  and  fatigue  crack  growth  laws  (see  bibliography  6-2) . 
Growth  of  disbonds  from  surface  flaws  or  from  edges  has  also  been 
studied.  In  all  such  studies,  strain  energy  release  rates  or 
stress  intensity  factors  are  calculated  by  analytical/numerical 
methods,  and  other  semiempirical  parameters  are  determined  from 
correlation  with  experimental  data.  In  most  of  these  experimental 
works,  the  growth  of  implanted  flaws  is  either  one-dimensional  or 
the  correlation  is  attempted  via  a  one-dimensional  model.  2-D 
finite  element  methods  (see  B. 2. 3. 6-12)  appear  to  be  an  attrac¬ 
tive  tool  for  numerical  calculations  because  of  their  versatil¬ 
ity  to  consider  various  boundary  conditions  like  free  edges  or 
surface  flaws.  The  methods,  however,  can  be  expensive,  especially 
when  their  accuracy  is  investigated  by  varying  mesh  sizes  and  ele¬ 
ment  shapes.  No  detailed  investigation  in  this  direction  dealing 
with  laminated  composites  is  reported  in  literature,  although  some 
studies  must  have  been  conducted  by  workers  who  have  used  such 
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techniques.  1-D  beam  theories  (B2. 1.4-5,  B2.2.2)  and  2-D  elas¬ 
ticity  solutions  (B2.2.1),  with  due  consideration  to  laminate 
stacking  sequence,  which  have  been  used  for  isolated  disbonds 
only,  are  comparitively  inexpensive,  but  extension  of  such  meth¬ 
ods  to  consideration  of  disbonds  near  free  edges  or  surface  flaws 
is  complicated  and  has  not  yet  been  attempted.  It  should  be 
pointed  out  that  techniques  of  modeling  self-similar  disbond 
growth  are  well  advanced,  although  a  significant  amount  of  ex¬ 
perimental/analytical  correlation  studies  is  needed  before  they 
can  be  used  as  practical  design  tools,  since,  in  many  cases,  the 
disbond  growth  is  not  self  similar.  Use  of  3-D  finite  elements 
has  also  been  attempted  for  calculation  of  strain  energy  release 
rates  (see  B2.1.7),  but  very  few  studies  are  reported  in  the  lit¬ 
erature  which  deal  with  two-dimensional  growth  of  disbonds.  Ob¬ 
viously,  such  finite  element  methods  are  attractive  tools  for 
studying  the  stress  field  near  disbonds  with  2-D  planforms,  but 
useful  results  for  such  problems  have  also  been  obtained  by  the 
use  of  modified  laminated  plate  theories  (ref.  3) .  No  study  on 
the  use  of  other  3-D  elasticity  theories  and  methods  has  yet  been 
attempted,  but  further  work  in  this  direction  is  needed  to  deter¬ 
mine  the  accuracy  of  the  techniques  which  are  currently  in  use. 
Last,  but  not  least,  a  considerable  amount  of  experimental  work 
is  still  necessary  to  determine  the  various  material  constants 
as  well  as  the  statistical  variability  of  those  parameters  which 
are  used  in  modeling  growth  of  disbonds. 


NADC- 8004 8-60 


EXPERIMENTAL  PROGRAM 


The  experimental  program  consisted  of  the  following  tasks: 

Task  I.  Study  the  effects  of  changes  in  disbond  position 
through  the  thickness  of  laminated  beams  on 
quasi-static  propagation. 

Task  II.  Study  the  effects  of  sharpness  of  the  disbond 

tips  on  quasi-static  failure  of  thick  laminated 
beams . 

Task  III.  Growth  of  disbonds  with  circular  planform  under 

cyclic  transverse  shear  in  thick  laminated  plates. 

TASK  I.  -  EFFECTS  OF  DISBOND  LOCATION 

Tests  were  performed  on  three  types  of  laminates  of  AS-3501 
graphite/epoxy  material  with  seven  different  defect  locations 
through  the  thickness  of  the  samples.  All  tests  were  performed 
in  three  point  bending  at  room  temperature. 

Test  Method 


Testing  was  performed  in  an  Instron  load  frame  at  a  rate  of 
1.27  mm  (0.05  inches)  per  minute.  Prior  to  testing,  the  static 
test  machine  was  balanced  and  calibrated  and  care  was  taken  to 
assure  that  the  force  was  introduced  perpendicular  to,  and  in 
the  center  of,  each  of  the  25.4  mm  x  254  mm  x  64  ply  graphite/ 
epoxy  beams.  Load  was  applied  using  a  12.7  mm  dia.  circular  steel 
rod  and  a  12.7  mm  wide  copper  pad  for  uniform  load  introduction 
and  to  avoid  stress  concentrations  at  the  center  of  the  beam. 
Geometry  of  the  test  specimens  is  shown  in  figure  1.  All  im¬ 
planted  defects  were  located  in  the  region  of  compressive  flex¬ 
ural  stress.  Ultrasonic  C-scans  were  made  on  all  beams  prior 
to  testing  and  if  catastrophic  failure  did  not  occur  after  a 
significant  load  drop,  the  beams  were  rescanned  and  retested. 
Examples  of  these  scans  are  given  in  figures  2  and  3. 
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Fabrication  and  Sample  Geometry 

Samples  were  made  from  AS-3501-6  graphite/epoxy  prepreg  in 
three  laminate  configurations  which  will  be  called  types  A,  B  and 
C  as  follows:  Type  A  =  [08/+454/-454 ] 2s ,  Type  B  =  [0g/+454/-454/ 

06//+454//“454y/02-'  s '  Type  c  =  ^°4//-452'/+452//°4  ^  2s  *  Al1  samPles  were 
64  ply  and  were  fabricated  with  25.4  mm  x  25.4  mm  two  ply  teflon 

defects  implated  in  locations  described  in  table  1.  Defect  loca¬ 
tion  is  shown  in  figure  1  and  was  verified  by  the  ultrasonic  C-scan 
technique.  Samples  are  identified  by  the  following  coded  example: 


C  -  20/21  -  2 


Specimen  number  2 

Defect  location  between  plies  20  and  21 
■Laminate  type  C  =  [04/+452/-452/04]  2S 


Three  replicates  for  each  location  were  tested  which  are  identified 
by  the  specimen  number,  i.e.  the  last  digit  in  the  example  shown 
above.  All  samples  were  cured  using  the  standard  cure  cycle  speci¬ 
fied  for  AS-3501-6  composites.  The  details  of  the  cure  cycle  are 
given  below. 

1.  Prepreg  laminate  is  prepared  with  bleeder  cloths  in  the 
ratio  of  one  layer  of  fiberglass  bleeder  to  three  plies  of 
graphite/epoxy . 

2.  Cool  plate  is  sealed  and  vacuum  is  drawn  to  -30  in.  Hg. 

3.  Temperature  is  increased  to  260°F. 

4.  Dwell  one  hour. 

5.  Pressurize  to  85  psi  with  Nitrogen. 

6.  Increase  temperature  to  350°F. 

7.  Cure  at  temperature  and  pressure  for  two  hours. 

8.  Decrease  temperature  to  below  250°  F. 

9.  Depressurize  and  cool  to  room  temperature. 
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Data 

All  samples  were  tested  until  a  significant  load  drop  or  failure 
occurred.  These  data  are  given  for  the  three  laminate  configurations 
in  tables  2,  3  and  4.  In  each  table,  the  following  data  are  given: 
sample  I.D.,  width,  thickness,  failure  load,  defect  lengths  1  and  2 
(as  viewed  in  the  C-scans)  and  defect  spacing.  Laminate  configura¬ 
tions  B  and  C  indicated  only  one  load  drop  coincident  with  failure . 
Laminate  type  A,  however,  experienced  as  many  as  four  load  drops 
prior  to  failure  and  these  data  are  given  in  table  5.  After  each 
load  drop,  the  samples  were  rescanned  and  tested  until  the  next  load 
drop  or  failure,  as  discussed  previously. 

TASK  II.  -  EFFECTS  OF  SHARPNESS  OF  DISBOND  TIPS 

Test  Method 

The  test  method  for  task  II  is  the  same  as  that  for  task  I . 

25.4  mm  x  254  mm  x  64  ply  graphite/epoxy  beams  with  defects  im¬ 
planted  between  the  mid-plies  were  loaded  in  three  point  bending. 

The  test  section  was  152.4  mm  long  and  the  beams  were  loaded  at 
the  center  of  the  test  section.  The  samples  were  staticly  loaded 
at  a  rate  of  1.27  mm  per  minute  in  an  Instron  load  frame.  Maximum 
load  was  recorded  on  a  chart  recoder  and  later  transferred  to  data 
tables.  Twelve  blunt  crack  specimens  and  12  sharp  crack  specimens 
were  tested. 


Fabrication  and  Sample  Geometry 

Twelve  samples  from  the  test  program  conducted  under  contract 
N62269-79-C-0209 ,  which  were  76.2  mm  x  254  mm  by  64  plies,  were 
used  to  make  the  sharp  crack  specimens.  These  plates  had  all  been 
previously  fatigued  in  three  point  bending.  Cyclic  loading  is  known 
to  yield  sharp  crack  fronts.  These  beams  were  all  ultrasonically 
C-scanned  to  determine  the  sizes  and  locations  of  defects.  The 
76.2  mm  wide  samples  were  marked  and  25.4  wide  beams  were  cut 
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from  the  center  region  of  each  plate.  These  beams  were  then  re¬ 
scanned  and  defect  sizes,  as  well  as  width  and  thickness,  were 
measured.  These  scans  and  data  are  given  in  figure  4  and  table  6. 
Average  defect  size,  radius  of  the  crack  front  and  crack  spacing 
were  determined.  These  data  were  used  to  fabricate  the  blunt 
crack  specimens.  The  blunt  crack  specimens  were  made  from  AS-3501- 
6  prepreg  and  the  delaminations  were  implanted  between  the  center 
plies  of  the  [ (04/±452/+452/04) 5] ^  laminates.  The  sample  geometry 
is  shown  in  figure  5,  where  the  defect  size  and  location  is  indi¬ 
cated.  The  average  defect  length  was  made  to  be  53.98  mm  with  a 
crack  front  radius  of  19.05  mm.  The  defects  were  made  from  two 
ply  teflon  as  in  task  I.  Specimen  identification  numbers  are 
prefixed  by  "B"  for  blunt  as  in  B-4  and  by  "S"  for  sharp  as  in 
S-12. 


Data 


Data  are  presented  in  table  6  which  contains  sample  I.D., 
thickness,  width,  disbond  lengths  1  and  2,  maximum  load  and  the 
defect  location  at  which  failure  occurred.  Ultrasonic  C-scans 
of  the  blunt  and  sharp  crack  specimens  are  shown  in  figures  4 
and  6,  respectively. 

TASK  III.  GROWTH  OF  CIRCULAR  DISBONDS  UNDER  CYCLIC  LOADING 

Test  Method 


Effects  of  size  and  location  of  disbonds  placed  in  [04/±452/ 
+452/04]2g,  AS-3501-6  graphic  epoxy  laminates  on  flaw  growth  un¬ 
der  cyclic  transverse  shear  were  studied  in  this  task.  This  was 
a  continuation  of  the  studies  on  disbonds  with  2-D  planforms  re¬ 
ported  in  reference  3.  Circular  disbonds  of  25.4  mm,  31.75  mm 
and  38.1  mm  diameter  in  the  form  of  teflon  films  were  placed  in 
76.2  mm  wide  x  254  mm  long  specimens.  Two  similar  disbonds  were 
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placed  symmetrically  in  the  net  section  span  of  152.4  mm  with 
their  centers  101.6  mm  apart.  All  fatigue  specimens  were  tested 
in  three  point  bending  in  a  manner  similar  to  the  beam  specimens 
in  Task  I  and  II.  Some  of  the  disbonds  were  placed  in  the  mid¬ 
plane,  i.e.  between  plies  32/33  and  others  were  placed  at  Z=±0.25H 
(between  plies  16/17  or  48/49)  ,  H  being  the  laminate  thickness. 

It  should  be  pointed  out  that  Z=+0.25H  indicates  that  the  disbonds 
were  in  the  region  of  compressive  flexural  stress  and  Z=-0.25H 
means  that  they  were  in  the  region  of  tensile  flexural  stress. 

The  samples  with  25.4  mm  diameter  defects  were  first  tested 
at  an  S-level  of  0.5  based  on  a  static  failure  load  of  29,400 
Newtons  as  determined  from  tests  reported  in  reference  3.  It 
should  be  noted  that  those  static  failures  were  due  to  flexure 
and  not  due  to  transverse  shear.  The  specimens  with  defects  lo¬ 
cated  at  Z=+0.25H  showed  little  propagation  after  12  million  cy¬ 
cles.  S-level  was,  therefore,  increased  to  0.75  for  31.75  mm  and 
38.1  mm  defects  located  at  Z=+0.25H.  S-level  was  kept  at  0.5  for 
all  other  samples.  The  ratio  R  of  minimum  to  maximum  load  was 
always  kept  equal  to  0.1.  Details  of  all  specimens  tested  are 
given  in  table  7,  where  measured  thicknesses  and  widths  are  also 
given.  The  specimens  are  identified  according  to  the  scheme  given 
below. 


3-1R0-1;  3  refers  to  task  III 

1R  refers  to  1  inch  (25.4  mm)  dia.  round  defect 
0  refers  to  Z=0,  midplane 
-1  refers  to  replicate  1 
3-125R-25-2;  125R  refers  to  1.25  inch  (31.75  mm) 
dia.  round  defect 
-25  refers  to  Z=- . 25H 

-2  refers  to  replicate  2 

3-150R+25-3;  150R  refers  to  1.5  inch  (38.1  mm)  dia. 
round  defect 
+25  refers  to  Z=+.25H 

-3  refers  to  replicate  3 
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Data 


On  an  average  the  25.4  mm  diameter  defect  samples  were 
scanned  about  24  times  during  their  life,  while  the  same  fig¬ 
ure  for  31.75  and  38.1  mm  diameter  defect  samples  was  about 
14.  Representative  C-scans  for  some  of  the  samples  are  given 
in  figures  7-15.  Growth  of  disbonds  at  three  locations  on 
four  fronts  labeled  as  A,  B,  C  and  D  (see  fig.  16)  were  mea¬ 
sured  from  each  of  the  scans  and  the  results  were  tabulated. 
These  data  were  later  used  for  data  correlation  studies  which 
are  reported  in  the  next  section. 
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DATA  CORRELATION  STUDIES 


Results  from  tasks  I  and  II  were  utilized  to  calculate  the 
critical  values  of  strain  energy  release  rates  with  the  help  of 
the  analytical  methods  developed  in  reference  2.  These  methods 
yield  the  stress  intensity  factors  and  energy  release  rates  at 
the  tips  of  two  symmetrically  located  disbonds  in  a  laminated  com¬ 
posite  beam  loaded  in  three  point  bending.  The  layers  are  assumed 
to  be  in  a  state  of  plane  stress.  Elastic  properties  of  the  layers, 
nondimensionalized  by  a  factor  EQ=6.895  GPa  (106  psi)  which  were 
used  for  calculations  are  given  below. 


Material 

Layers 

cn/Eo 

C33/E0 

C13/E0 

C55/E0 

1 

0° 

18.24 

1.46 

0.409 

0.841 

2 

+45°  or  -45° 

2.16 

1.49 

0.294 

0.731 

3 

(+45/-45J - 
smeared 

2.88 

1.50 

0.130 

0.687 

4 

Epoxy  Matrix 

0.605 

0.605 

0.228 

0.188 

TASK  I.  -  EFFECT  OF  DISBOND  LOCATION 

Since  the  +45°  and  -45°  layers  are  lumped  separately  in  laminate 
types  A  and  B,  properties  listed  under  material  no.  2  were  used  for 
both  of  these  layers.  For  consistency,  the  same  properties  were  also 
used  for  laminate  C.  Calculations  indicate  that  if  the  properties 
of  a  smeared  material  (no.  3)  are  used,  the  results  differ  by  less 
than  five  percent. 

The  computer  code  DELAM,  in  its  present  form,  can  solve  the 
problem  of  disbonds  between  two  layers  of  the  same  material.  In 
a  majority  of  test  specimens,  this  was  the  case  and  60  nonzero  terms 
of  the  Fourier  Series  were  enough  for  the  calculation  of  the  regular 
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kernels  of  the  integral  equations.  Forty  collocation  points  yielded 
reasonable  accuracy  in  calculation  of  the  unknown  functions,  inten¬ 
sity  factors  and  energy  release  rates.  In  a  few  cases,  however,  the 
disbonds  were  between  a  0°  and  a  +45°  (or  a  -45°)  layer.  Although, 
the  code  can  be  modif ed  to  handle  these  cases ,  this  was  beyond  the 
scope  of  the  current  program.  Therefore,  calculations  were  performed 
for  disbonds  located  between  two  epoxy  matrix  layers  (material  no.  4) 
of  relatively  small  thicknesses  which  are  sandwiched  between  0°  and 
45°  materials.  In  practice,  there  is  always  a  matrix  rich  region 
between  two  laminae.  Results  were  obtained  for  three  gradually  de¬ 
creasing  values  of  thickness  t  of  the  matrix  layers  (up  to  t^/t  = 
.000734,  l  being  the  span  of  the  beam)  and  tho  enery  release  rates 
in  the  limiting  case  of  zero  thickness  were  obtained  by  extrapolation. 
A  three  point  extrapolation  formula  worked  well  for  calculation  of 
Gi;r  as  it  appeared  from  results  of  the  cases  of  disbonds  between  two 
similar  layers.  For  calculation  of  Gj,  however,  the  extrapolation 
procedure  was  not  very  satisfactory,  especially  when  the  calculated 
values  of  Gj  were  of  the  order  of  Gjj.  This  phenomenon  is  expected 
since  the  nature  of  the  stress  singularity  changes  from  the  inverse 
square  root  type  to  the  type  -1/2  ±  ie,  in  the  limiting  case.  More¬ 
over,  as  the  thickness  of  matrix  layers  is  gradually  reduced,  it  is 
necessary  to  take  increasing  number  of  terms  in  the  Fourier  series 
for  calculation  of  the  regular  kernels.  Three  hundred  nonzero  terms 
were  considered  for  the  smallest  thickness .  It  also  appeared  that 
further  reduction  in  the  thickness  may  necessitate  the  use  of  more 
collocation  points  (>40)  .  For  cases  where  Gj  was  of  the  order  Gjj 
the  values  of  G^.  showed  a  little  oscillatory  tendency  about  an  aver¬ 
age  value.  In  such  cases,  the  limiting  value  of  Gj  was  taken  as  the 

average  of  the  values  calculated  for  three  values  of  t  ,  i.e.  t  /£  = 

m  m 

.000734,  .001468  and  .002202  was  utilized.  Considering  the  scatter 
in  experimental  data,  these  results  appeared  satisfactory  for  the 
purpose  of  data  correlation  studies  reported  herein. 
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Calculated  values  of  GjEqJ l/p2,  GjjEqVP2,  ( £  =  span  of  the  beam 
and  P  =  applied  load  per  unit  width)  as  well  as  the  total  energy 
release  rates  are  shown  in  tables  8,  9  and  10  for  laminate  types  A, 

B  and  C,  respectively.  Use  of  these  values  and  the  measured  load 
carrying  capacities  of  the  specimens  yielded  G ^  the  critical  values 
of  energy  release  rates,  which  are  also  shown  in  the  tables.  It 
should  be  noted  that  for  laminate  types  B  and  C  the  failure  loads 
were  also  the  maximum  loads.  Some  specimens  of  type  A  laminates, 
however,  showed  some  load  drops  before  final  failure  occurred  as  can 
be  seen  in  table  5.  In  some  cases,  i.e.,  for  some  of  the  specimens  of 
the  Series  A-16/17,  A-28/29  and  A-32/33,  the  first  load  drop  occurred 
at  very  low  loads.  The  reasons  for  such  load  drops  before  failure 
are  not  very  clear.  It  should  be  pointed  out,  however,  that  the 
present  analysis  does  not  consider  the  effect  of  the  tearing  mode 
(GII;I)  ,  which  can  have  significant  influence,  especially  in  cases 
where  the  +45°  and  -45°  layers  are  lumped  together.  Also,  shifting 
of  the  disbond  location  across  the  thickness  of  the  specimen  may 
yield  such  load  drops.  These  effects  can  not  be  assessed  at  this 
time  but  should  be  considered  in  future  studies.  For  the  current 
study,  the  value  of  the  maximum  load  was  used  for  calculations. 

Since  the  disbond  geometry  might  have  changed  before  the  maximum 
load  was  attained,  the  results  for  laminate  type  A  should  be  inter¬ 
preted  with  some  caution. 

All  calculations  were  performed  for  a  laminate  thickness  of 
8.95  mm.  and  a  beam  width  of  25.4  mm.  Experimental  values  of  loads 
were  normalized  with  due  consideration  to  specimen  thickness  and 
width  and  these  values  are  shown  in  tables  8-10.  Considering  the 
3catte  •  of  experimental  data  which  exists  in  results  for  three  test 
specimens  of  the  same  type  and  values  of  Gjj  reported  in  table  11 
(results  of  task  II  discussed  in  the  next  subsection) ,  it  appears 
that  the  critical  energy  release  rate  does  not  vary  to  a  great  ex¬ 
tent  with  its  location,  i.e.  the  orientation  of  the  layers  between 
which  the  disbond  is  located.  This  is  expected  if  one  considers  the 
fact  that  there  is  always  a  matrix  rich  region  between  two  laminae 
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and  contends  that  Gc  is  a  characteristic  material  constant  for  the 
matrix  material,  which  may,  of  course,  depend  on  its  type  and  proces¬ 
sing  variables. 

In  a  few  cases,  A-12/13,  A-16/17,  B-6/7,  B-12/13,  C-8/9  and 
C-12/13,  calculated  values  of  GT  is  of  the  order  of  GTI  and  the 
values  of  G  ,  obtained  in  these  cases  are  consistently  on  the  low 
side.  This  suggests  that  a  quadratic  interaction  criterion  of  the 
type 


II.  2 
.c 
’II 


*  1 


(1) 


instead  of  the  one  used  here,  Gj  +  G -  Gc,  may  yield  a  better  fit 

to  the  experimental  data.  If  G^  is  taken  equal  to  1040  N/m  (as  de¬ 

termined  in  task  II,  see  table  11) ,  calculations  indicate  that  G^ 
should  lie  somewhere  between  200  to  400  N/m.  No  detailed  examination 
of  such  a  criterion  is  attempted  here,  since  the  scatter  in  pres¬ 
ent  data  on  G^  does  not  yield  a  realistic  assessment  of  G^ .  Alter¬ 
native  experiments  with  Gx  dominating  the  response,  should  be  designed 
for  this  purpose. 

Values  of  Gc  obtained  from  the  specimens  of  the  series  B-26/27 
are  a  bit  on  the  high  side,  reasons  for  which  are  not  quite  clear. 

In  this  case,  the  disbond  is  located  between  +45  and  -45  and  effects 

of  mode  III,  which  is  not  considered  in  the  analysis,  may  have  some¬ 

thing  to  do  with  these  results.  It  is  noted,  however,  that  other 
specimens  with  disbonds  between  +45  and  -45  layers,  i.e.  B-12/13, 
A-12/13  and  A-28/29  yield  values  of  GC  which  are  on  the  low  side. 

It  should  be  pointed  out  that  in  cases  when  Gj  +  0,  the  stress 
intensity  factor  Kj  at  one  of  the  tips  is  positive,  whereas  at  the 
other  tip,  it  is  negative,  indicating  existence  of  compressive 
stresses  near  that  tip.  Therefore,  it  is  quite  possible  that  the 
two  surfaces  of  the  disbond  may  come  in  contact  near  such  a  tip 
causing  frictional  resistance.  The  effect  of  this  phenomenon  cannot 
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be  assessed  here,  but  obviously,  it  may  have  some  influence  on 
the 'results  presented.  Tables  8-10  show  the  location  of  the  tip 
where  Kj  was  found  to  be  positive.  Location  c  indicates  the  tip 
nearer  to  the  supports,  and  d  identifies  the  tip  near  the  applied 
load. 

To  examine  the  usefulness  of  the  analytical  methodology  in 
assessing  criticality  of  disbonds  placed  at  different  locations, 
failure  loads  P-  for  the  three  types  of  laminates  predicted  by 

the  criterion  G  ■  G_  +  GTT  ■  G  =  1000  N/m  are  plotted  in  figures 

^  2  c 

17-19..  It  may  be  noted  that  the  criterion  yields  Pf  =  G  EQi/a 

where  the  factor  a  =  GEoi/P2,  listed  in  tables  8-10.  The  figures 
show  reasonable  correlation  with  experimental  data  which  have  been 
normalized  for  2.54  cm  wide  beams  with  2.54  cm  size  defects.  In 
general,  the  differences  between  predicted  and  experimental  values 
are  higher  when  contributions  of  G ^  to  the  total  energy  release 
rates  are  relatively  high,  indicating  that  the  quadratic  inter¬ 
action  criterion  discussed  earlier  or  even  the  assumption  of  lin¬ 
ear  interaction  with  different  coefficients  associated  with  Gj 
and  GI;[  may  yield  better  correlation  with  experimental  data. 

TASK  II.  -  EFFECT  OF  SHARPNESS  OF  DISBOND  TIPS 

Smeared  (+45 /—  45)  material  properties  (material  no.  3)  were 

2  2 

used  for  calculation  of  G^EqJI/P  =  GEQ£/P  (GI=0)  for  disbonds 
in  laminates  with  a  thickness  8.95  mm  and  a  width  of  25.4  mm  of 
the  type  tested  in  task  II.  Normalized  values  of  failure  loads, 
with  due  consideration  to  width  and  thickness  of  the  specimens 
were  used  for  determination  of  Results  are  given  in  table 

11.  Results  show  some  scatter  and  about  10%  difference  in  aver- 
age  Gjj  for  blunt  and  sharp  tips.  Considering  the  scatter,  this 
difference  is  not  significant  and  it  appears  that  the  nature  of 
the  tip  does  not  have  any  dominant  influence  on  catastrophic  dis¬ 
bond  growth  under  static  shear. 


NADC-80048-60 


TASK  III.  -  PROPAGATION  OF  DISBOND  IN  FATIGUE 


Strain  energy  release  rates  along  the  periphery  of  ellip¬ 
tic  disbonds  of  various  sizes  with  semi-axes  and  L 2  located 
at  the  midsurface  of  [0^/±452/+452/0^ ] 2 s  (which  can  also  be  writ¬ 
ten  as  (  (0^/±452/^452/0^) g]  )  AS-3501-6  graphite  epoxy  laminates 
were  reported  in  reference  3.  These  results  were  also  utilized 
in  the  present  study.  Properties  of  the  unidirectionally  rein¬ 
forced  composite  (0°  layer)  and  stiffness  of  the  (04/±452/t452/ 
04)g  laminate  are  reproduced  here  in  table  12.  As  indicated  in 
reference  3,  D^g  and  D2g  were  neglected  for  the  calculation  and 
the  laminate  was  considered  as  orthotropic.  Since  the  disbond 
is  located  in  the  mid-plane,  the  energy  release  rate  G^.  in  open¬ 
ing  mode  is  identically  equal  to  zero  in  this  case.  In  the  pres¬ 
ent  study,  disbonds  were  also  placed  at  Z=±0.25H,  H  being  the 
laminate  thickness.  For  such  location,  Gj  is  nonzero  and  the 
variations  of  GII#  Gm  and  Gi  alon9  the  disbond  periphery  are 
shown  in  figures  20-22.  The  calculations  were  performed  with 
laminate  properties  given  in  table  13,  for  self  equilibrating 
uniform  shear  stress  in  x-direction  on  the  disbond  surfaces. 


The  two  laminates  were  considered  orthotropic,  i.e.  Ai6=A26=D16= 
°26=®  i°r  purpose  of  calculations.  As  in  the  case  of  disbond 
in  the  mid-plane,  G.^  first  increases  with  and  then  starts 
decreasing  provided  L2  is  kept  fixed.  Similar  behavior  is  also 
observed  for  Gj.  Figure  23  shows  the  variation  of  G^  and  Gj 
at  x  »  Lp  y  ~  0  with  L^/L2  for  fixed  values  of  L^.  It  should 
be  noted  that  in  the  limiting  case  L^/L^  +  0,  the  problem  re¬ 
duces  to  that  of  generalized  plane  strain.  Results  from  the 
2-D  elasticity  solution  (ref.  2)  for  this  problem  are  also  shown 


in  the  figure.  Although  the  results  for  this  limiting  case  can¬ 
not  be  obtained  by  the  methods  used  here,  the  results  indicate 
that  G^  is  underestimated,  as  was  observed  in  reference  3  for 
disbonds  in  the  mid-plane  of  the  laminates.  It  is  clear  from 
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the  results  given  here  that  Gj  is  overestimated  by  the  laminated 
plate  theory  assumptions.  Total  energy  release  rate  (Gj  +  G  ) 
is,  however,  underestimated. 

Growth  of  disbonds  under  cyclic  loading  observed  experimen¬ 
tally  were  measured  and  tabulated  as  described  in  the  previous 
section.  Growth  in  y  direction  (in  the  direction  of  the  width 
of  the  plate  specimens)  was  not  significant.  Growth  rates  in 
x-direction  were  measured  at  6  points  (y=0  and  y=±9.53  mm).  The 
distance  1^  of  the  two  disbond  fronts  at  y=0  from  the  center  of 
each  of  the  two  flaws  were  plotted  against  N,  the  number  of  cy¬ 
cles  for  each  specimen.  Some  representative  data  are  shown  in 
figures  24-28.  Growth  rates  AL^/AN  were  then  computed  at  various 
points  along  the  L1  vs.  N  plots.  Values  of  the  current  total 
length  of  the  disbond  in  x-direction  were  then  used  as  the  dimen¬ 
sion  of  the  major  axis  of  an  approximated  elliptic  disbond  to 
calculate  A(/G^j,  the  range  of  variation  of  /G^ .  Values  of 
dL1/dN  were  then  plotted  against  A(/G^)  on  log- log  graph  papers. 

Results  for  disbonds  in  mid-plane  are  shown  in  figure  29. 
Least  square  fits  were  attempted  to  obtain  the  empirical  con¬ 
stants  and  n  in  the  relationship  of  the  form: 


dL 


N  C1  (/GII)  “  A0  } 


n 


(2) 


A  series  of  values  of  the  threshold  range  An  (/gtt) ,  below  which 

u  J"L  1/2 

flaw  growth  is  not  expected  varying  between  7  to  11  (N/m)  were 
considered.  The  constants  for  two  curves  which  appear  to  yield  good 


are  given 

below  in 

Newton  meter  units. 

(N/m) 1/2 

C^(N/m  units) 

n 

00 

• 

o 

0.1456  x  10“10 

4.032 

10.0 

0.2782  x  10"9 

3.095 
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The  curve  corresponding  to  equation  (2)  with  the  second  set  of 
parameters  is  shown  in  figure  29.  These  parameters  differ  a 
little  bit  from  those  obtained  in  reference  3  for  disbonds  in  a 
similar  location,  i.e.  mid-plane.  The  difference  may  be  due  to 
difference  in  mean  stress  or  S  values.  S  values  for  the  speci¬ 
men  tested  under  the  present  program  were  0.5,  whereas  for  those 
considered  in  refernce  3,  S  varied  between  0.4  to  0.6.  It  should 
also  be  pointed  out  that  for  the  maximum  load,  values  of  Gjj  at¬ 
tained  is  of  the  order  of  and,  therefore,  fitting  of  curve  of 
the  type  (2)  may  not  be  appropriate  for  such  high  values  of  GI3;. 
Curve  fits  to  consider  differing  S  values  and  high  values  of  G^ 
are  currently  being  attempted  in  an  ongoing  follow-on  effort. 

Results  for  disbonds  located  at  Z=±0.25H  are  shown  in  figure 
30  along  with  the  curve  fit  to  data  for  disbonds  at  Z=0,  described 
above.  It  is  noted  that  these  test  data  are  for  two  values  of  S, 
i.e.  S=0.5  and  S=0.75.  Since  the  static  failure  load  was  deter¬ 
mined  from  tests  which  indicated  flexural  failure  and  the  shear 
stresses  at  mid-plane  and  at  Z=±0.25H  for  a  fixed  load  are  not 
the  same,  the  actual  mean  stresses  for  specimens  with  disbonds 
in  mid-plane  and  for  those  with  defects  at  Z=±0.25H  could  be 
different.  Also,  for  defects  located  Z=±0.25H,  the  problem  is 
more  complicated  because  of  the  exitence  of  Gj,  which  may  increase 
the  flaw  growth  rate.  The  values  of  A(i/G^)  for  these  specimens 
were  about  half  of  those  of  A(/S^).  An  equation  to  consider 
simultaneous  influence  of  A(/G^)  and  A(/G~),  may  be  postulated 
but  such  an  attempt  is  not  made  here.  It  appears,  however,  that 
all  test  data  for  S=0.5  (all  of  these  points  be  on  the  left  hand 
side  of  fitted  curve)  show  a  slightly  higher  growth  rate  than 
disbonds  at  mid-plane.  The  data  points  for  S=0.75  on  the  other 
hand  fall  close  to  a  vertical  line  indicating  a  kind  of  unstable 
growth  pattern,  possibly  because  of  the  fact  the  maximum  values 
of  G  attained  in  these  tests  are  close  to  that  needed  for  cata¬ 
strophic  growth  under  quasi-static  loading. 
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Before  closing  this  discussion,  it  should  be  pointed  out 
that  growth  of  a  disbond  under  pure  mode  II  loading  should  be 
the  same  at  x^iL^  as  predicted  by  the  empirical  growth  law  and 
the  stress  analysis  model  used  here.  The  results  for  beams 
with  mid-plane  disbonds,  however,  show  some  differences  between 
the  growths  at  the  four  fronts,  possibly  because  of  the  presence 
of  inhomogeneities  and/or  the  approximations  used  in  the  analy¬ 
sis.  This  is  one  of  the  reasons  for  the  scatter  noticed  in  data 
correlation.  Presence  of  mixed  mode  conditions  for  dibonds  at 
Z=+0.25H  may  introduce  larger  differences  because  of  the  exis¬ 
tence  of  tensile  and  compressive  normal  stresses  at  the  two 
fronts.  As  mentioned  before,  better  curve  fits  to  fatigue  data 
are  currently  being  attempted  in  an  ongoing  program. 


NADC-8004  8-60 


DISCUSSIONS  AND  CONCLUSIONS 


In  the  present  work  the  state  of  the  art  of  assessing  de¬ 
fect  criticality  in  laminated  composites  has  been  summarized  and 
the  significance  of  various  ongoing  programs  in  this  area  have 
been  assessed.  In  addition,  a  number  of  tests  have  been  conducted 
and  data  obtained  from  these  tests  have  been  analyzed  to  test  the 
range  of  validity  of  the  criticality  criteria  developed  in  re¬ 
lated  past  programs.  From  the  results  reported  herein  the  fol¬ 
lowing  conclusions  and  comments  can  be  made. 

1.  A  linear  elastic  fracture  mechanics  approach  and 
methods  of  stress  analysis  based  on  2-D  elasticity  theory 
can  be  used  to  assess  criticality  of  disbonds  in  laminated 
beam  type  structures  under  transverse  shear,  irrespective 

of  disbond  size  and  location.  Validity  of  the  proposed  meth¬ 
odology  has  been  demonstrated  for  cases  where  mode  II 
energy  release  rate  is  the  dominant  factor  influencing 
failure.  Further  work  appears  necessary  to  assess  the 
mode  III  energy  release  rate  as  well  as  the  interaction  of 
all  three  modes;  i.e.  I,  II  and  III,  in  causing  failure. 

2.  Sharpness  of  disbond  tips  introduced  by  fatigue 
loading  does  not  appear  to  have  a  significant  influence 

c 

on  Gjj,  as  seen  from  comparison  with  test  data  from  dis¬ 
bonds  with  blunt  tips,  i.e.  implanted  teflon  defects. 

However,  because  of  the  scatter  in  experimental  results, 
more  tests  are  needed  before  a  final  conclusion  is  drawn. 

3.  Methods  of  stress  analysis  based  on  laminated 
plate  theory  yield  reasonably  accurate  values  of  energy 
release  rates  along  the  periphery  of  disbonds  in  laminated 
plate  type  structures  under  transverse  shear.  Usefulness 
of  semi-empirical  crack  growth  laws  in  predicting  the 
growth  of  such  flaws,  for  cases  where  mode  II  energy  re¬ 
lease  rate  is  dominant,  has  been  clearly  illustrated. 

Further  work  for  evaluating  effects  of  varying  mean  stress 
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as  well  as  influence  of  mode  I  and  mode  II  energy  release 
rate  on  flaw  growth  are  needed  to  put  the  NDE  methodology 
for  assessing  criticality  of  disbonds  under  fatigue  load¬ 
ing  on  a  sound  basis. 
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Table 

7.  Details 
Tested 

of  Samples  with  Circular  Disbonds 
in  Fatigue  -  Task  III 

Sample 

S-Level 

Max.  No.  of 

Cycles,  103 

Thickness , 
mm. 

Width 
mm . 

3-1R0-1 

.5 

380 

9.68 

76.0 

-2 

520 

9.14 

74.5 

-4 

210 

8.84 

75.8 

3-1R+25-1 

.5 

2,952 

9.53 

76.0 

-2 

11,345 

9.19 

76.4 

-3 

750 

9.50 

75.9 

3-1R-25-1 

.5 

290 

8.31 

75.9 

-2 

680 

8.66 

76.3 

_ -3 

940 

8.74 

76.3 

3-125R0-1 

.5 

36 

8.79 

75.6 

-2 

28 

8.64 

76.0 

-3 

16 

8.31 

75.7 

3-125R+25-1 

.75 

27.5 

8.59 

75.9 

-2 

5.55 

8.28 

75.8 

3-125R-25-1 

.5 

32 

9.75 

76.0 

-2 

32 

9.70 

76.4 

-3 

30.5 

8.92 

76.3 

3-150R0-1 

.5 

7 

7.59 

75.3 

-2 

14.77 

8.26 

74.9 

-3 

15 

8.28 

75.8 

3-150R+25-1 

.75 

4.61 

8.61 

75.1 

r 

p 

B 
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Table  7.  Details  of  Samples  with  Circular  Disbonds  (continued) 

Tested  in  Fatigue  -  Task  III 


Sample 

S-Level 

Max.  No.  of 
Cycles,  103 

Thickness, 

mm. 

Width 

mm. 

3-150R-25-1 

.5 

18 

8.31 

74.0 

-2 

18 

8.36 

75.6 

-3 

16 

9.68 

76.4 

Table  8.  Data  Correlation  for  Type  A  Laminates  -  Task 
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Table  8.  Data  Correlation  for  Type  A  Laminates  -  Task  I  (continued) 
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Table  9.  Data  Correlation  for  Type  B  Laminates  -  Task  I  (continued) 
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Table  10.  Data  Correlation  for  Type  C  Laminates  -  Task 
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Table  10.  Data  Correlation  for  Type  C  Laminates  -  Task  I  {continued) 
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Table  11.  Calculated  Values  of  Critical  Energy  Release  Rates 

for  Blunt  and  Sharp  Tips 


NADC-80048-60 
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Table  12.  Elastic  Properties  of  0°_Layer  and 
Stiffnesses  of  (0 ./±455/+45_/0 . ) 
Laminate  *  1 


Transversely  Isotropic  0°  Layer 


Ea  =  125  GPa 

Et  =  10  GPa 

v_  =  0.28 
A 

G.  =  5.8  GPa 
A 


(0 ./±45_/+45_/0.)  Laminate 

4 _ 2 _ 2  4  S _ 

A11  =  4,06  X  N/Jn* 

A12  =  0,79  x  1C)8  N/m- 
A22  =  1.24  x  108  N/m. 

Agg  =  0.95  x  108  N/m. 
A16  =  A26  °’° 

Ba6  =  °-° 
K55  =  2.12  x  107  N/m. 


D11  =  877  Nm* 

D12  =  131  Nm. 

D22  =  218  Nm. 

Dgg  =  160  Nm. 

D^g  =  ^2g  =  3.25  Nm. 

(a, 6  =  1,2,6) 

K44  =  1.75  x  107  N/m 


Ply  Thickness  =  0.1519  mm 
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Table  13.  Stiffnesses  of  (0 ./+45_/-45_/0 . )  and 
(047+452/^452/04)2  Laminates 


(°4/t452/?452/04) 

AX1  =  2.03  x  108  N/m 
A12  =  0.394  x  108  N/m 
A22  =  0.623  x  108  N/m 

O 

Agg  =  0.466  x  10  N/m 
A16  =  A26  =  0 

Ba$  “  0 

K  =  1.06  x  107  N/m 
s  s 


D11  =  138  Nm 

D12  =  7,4  ^ 

D22  =  16.7  Nm 

Dgg  =  11.0  Nm 
D16  =  D26  =  1*63  Nm 
(a, 8  =  1,2,6) 

K44  =  0.875  x  107  N/m 


(04/±452/+452/04)3 

A11  =  6.09  x  108  N/m 
A^2  =  1.18  x  108  N/m 
A22  =  1.87  x  108  N/m 
A66  =  1.40  x  108  N/m 

A16  "  A26  “  0 


2817  Nm 
488  Nm 
787  Nm 
585  Nm 
D26  =  1.63  Nm 


B  Q  *  0 

aS 


(a, S  =  1,2,6) 


K55  =  3.18  x  107  N/m 


K 


44 


2.63  x  10  N/m 


Ply  thickness  ■  0.1519  mm 


Point  Bend  Tests 
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Representative  C-scans  After  Load  Drop  and  Prior  to  Complete  Failure 


Figure  5.  Blunt  Crack  Sample  Geometry 
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Figure  8.  Flaw  Growth  in  Sample  3-125RO- 


Figure  9.  Flaw  Growth  in  Sample  3-150RO- 


cles 


Figure  12.  Flaw  Growth  in  Sample  3-125R+25-1,  S  =  0.75 
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Figure  15.  Flaw  Growth  in  Sample  3-150R-25- 


N ADC-80048-60 


Figure  16.  Measurement  Locations  for  Four  Crack  Fronts 


+  Predicted 
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Figure  18.  Failure  Loads  for  Disbonds  of  2.54  cm.  Length  in  2.54  cm.  Wide  Beams-Type  B  Laminates 


.895  GPa 


Figure  20  .  Variation  of  G  with  Dimensions  for  Flaws  at  Z  =  +0.25H 


Figure  21.  Variation  of  GII3;  with  Dimensions  for  Flaws  at  Z  =  ±0.2511 


Figure  22.  Variation  of  GT  with  Dimensions  of  Flaw  at  Z  =  ±0.25H 
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Figure  27. 


Figure  28.  L.  vs.  N,  Sample  3-150R-25- 
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